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Abstract 

The method of investigation and numerical solution of 
synthesis problems of antennas with flat radiating aperture 
according to the prescribed requirements to amplitude 
directivity pattern (DP) is proposed in this paper. Freedom of 
choice of the phase DP is used as an additional degree of 
freedom to improve the quality of approximation of the 
amplitude of the synthesized DP to the given one. It is shown 
that the non-uniqueness and branching of solutions, 
dependent on the size of two physical parameters that 
describe the size of antenna aperture, and the properties of 
given amplitude DP is characteristic features of this class of 
problems. The method of finding the branching lines of 
solutions of basic equations of synthesis, which is based on 
implicit function method, is proposed. This method allows to 
reduce a non-linear two-parameter spectral problem to 
solution of the corresponding Cauchy problem for 
differential equations of first order. On this basis the 
methodology of localization of existing solutions is 
developed. Effective computational algorithms for finding 
the optimal solutions of synthesis problems are constructed, 
the numerical examples are shown. Analysis of the 
effectiveness of different types of existing solutions for 
synthesis of several types of given amplitude DP is presented. 
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Introduction 

As is known [1], in many practical applications at the 
design stage of antennas the requirements are only to 
the energy characteristics of antenna (amplitude DP or 
DP by the power). Here freedom of choice of the phase 
DP can be used to improve the quality of 



approximation of synthesized DP to the given one. 
The synthesis problems of antennas according to the 
prescribed amplitude DP or DP by power belong to 
the class of problems with incomplete input 
information. Nonuniqueness and branching (or 
bifurcation) of existing solutions is characteristic 
feature of this class of problems. 

The problem of nonuniqueness and branching of 
existing solutions, determination of their general 
structure and main properties is investigated only 
partly for the synthesis problems of linear antennas 
and antenna arrays [2, 3]. It is shown that the quantity 
of existing solutions and their properties depend on 
the physical parameter c characterizing the length of 
antenna and solid angle, in which the prescribed 
amplitude DP is given, and properties of this DP. It is 
found that for arbitrary (bounded) nonnegative 
values of the parameter c there exist several types of 
solutions, belonging to the class of in-phase DP (called 
as primary). These solutions are optimal only at small 
values c . With the increase of the values c there exist 
such values c i the points of branching in which the 
more effective complex solutions branch-off from the 
primary (real) solutions. 

Note, the results obtained in [2, 3] may not be 
applicable to the synthesis problems of antenna with 
flat radiating aperture, since in these problems the 
solutions depend on two physical parameters and in 
contrast to the branching points (specific to the 
problem of synthesis of linear antennas and antenna 
arrays) there exist the branching lines of solutions. 
Moreover the problem of finding them is 
insufficiently studied nonlinear two-parameter 
spectral problem. The existence of connected 
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components of the spectrum, which in the case of real 
parameters are spectral lines, is an essential feature of 
nonlinear two-parameter problems. 

In mathematical aspect in variational statements the 
synthesis problems of antenna with flat radiating 
aperture are reduced to study two-dimensional 
nonlinear integral equations of the Hammerstein type. 
Accordingly to the general branching theory of 
solutions corresponding linear homogeneous integral 
equation with a nonlinear occurrence of two spectral 
parameters in the kernel of operator to find the 
branching lines is obtained. 

A numerical method for solving the nonlinear two- 
parameter spectral problems that allows to find the 
branching lines of solutions is proposed and justified in 
the work. Methods of the general branching theory of 
solutions allow to find the branching-off solutions. 
Jointly this gives opportunity to locate the existing 
solutions of the basic equations of synthesis and build 
effective calculation algorithms for finding and 
analysis of optimal solutions. 

Formulation of problem. Basic equations of 
synthesis 

According to [2] we shall formulate the synthesis 
problem of antennas with a flat radiating aperture G 
putting, that radiator is acoustic flat antenna placed in 
the plane XOY , and its DP is described by the formula 
[4] 

/(9,cp) = AU = JJ£/(jc, y)e *(^°^y™^) dxdy _ (1) 

G 

Here U(x,y) is the distribution of oscillating 
velocity on a 

flat surface of antenna, k = 2n/X is a wave number, X 
and is the wave length. 

Introduce the generalized angular coordinates 
s l = (sin8cos(p)/siny[ , s 2 = (sin8sin(p)/siny 2 (2) 

and the corresponding nondimensional parameters: 
q = ka x sin y 1 , c 2 = ka 2 sin y 2 , (3) 

setting that antenna (aperture G ) (Fig. 1) is in some 
rectangle G P = i\x\<a l , |y|<a 2 }, the geometric center 

of which coincides with the beginning of coordinates 
and the sides are parallel to the axes OX , OY . Put that 
amplitude DP F(Q,(p) is given also in the domain Q 

belonging to some 'rectangle' Q P = {\s { \ < b v |s 2 |<fr 2 } 



whose sides are parallel to the coordinate axes of the 
generalized coordinate system (2). The intervals of 
change of the angle 9 in which amplitude DP 

i Z 




FIG. 1 TO STATEMENT OF THE SYNTHESIS PROBLEM OF FLAT 
APERTURE 

F(s l ,s 2 ) is given at(p = cp = ji/2 we denote as yj 
and y 2 , respectively. The parameters c x and c 2 
characterize the dimensions (in wave lengths) of 
antenna (aperture) and domain (solid angle) Q , 
where amplitude DP is given. 

Using new variables, (1) takes the form 

/(«!, s 2 ) =AUs jju(x, y)e i{wc ^ ) dxdy , (4) 

G 

which we shall consider as the action of linear integral 
operator A from complex space of square integrable 
functions H v = (G) , describing the distribution of 
the oscillating velocity, into complex space of square 
integrable functions H f = l^iM. 2 ) . The properties of 

the operator A are defined by the type and geometry 
of the radiating system. 

We define scalar products and norms in the spaces 
H„ and H f as follows: 



{U l ,U 2 ) H =^-\\ U m ( X ,y)-U (2 \x,y)dxdy , 



R 2 

=(/./)£• (6) 
For representation (4) the Percival equality is valid [5]: 



30 



Advances in Applied Acoustics Vol. 1 Iss. 1, November 2012 



\ AU t=\\ u \L 



(7) 



From here follows the property of isometricity of the 
operator A . On the basis of introduced scalar products 
and relation (AU, f) Hf = (U, A* f) Hu we obtain, 

necessary later on, the expression for conjugate 
operator 



A*f = jjf(s l ,s 2 )e- i ^ +c ^ ) ds 1 ds 2 



(8) 



Consider the synthesis problem of given amplitude DP 
F(s l ,s 2 ) in the domain Q as a minimization problem 
of the functional 

°F ( U ) = jjj [ F ( s V s 2)-\f( s U s 2)\] 2 |*i*2 + 

+ j| |/(Ji,J 2 )| 2 *1*2 (9) 

in the space H v = (G) . Here the first summand 
describes the mean-square deviation of the modulus 
given and synthesized DP in the domain Q . The 
second summand imposes constraints on the level of 
energy emitted outside given solid angle (the 
domain Q ). 

Differentiating functional a F (U) by Hato and using 
the necessary condition of minimum functional [6] 
gmd<5 F (U) = 0, we write equation with respect to the 
function of distribution of the oscillating velocity in an 
aperture: 

U(x, y)= C ' C2 \\ F(s l ,s 2 )exp(-i(c l xs l +c 2 ys 2 ))x 
{2nf 3 l 



<exp 



i arg JJl/(je\ y'W^'^'^dx'dy' 



ds\ ds 2 . (10) 



The operator form reads: 

U = ^\A\F -Qx V (i3XgAU)) . 
An 



(11) 



Acting on both sides of this equation by the operator 
A and taking into account (4) we obtain the equation 
with respect to the synthesized DP: 

/(5 1 ,. 2 ) = B/-ffFU,4)^ 1 ,5 2 ,5|,4;c 1 ,c 2 ) e iaig/(s '' s i ) ^4 7 (12) 



K(Q,Q',c) = -^\\\exp[i(c lX ( S [ - Sl ) + c 2 y(s' 2 - j 2 ))] dxdy (13) 
(2ji) J ^ 

is a kernel essentially dependent on the form of the 
domain G . 

In the case of symmetric domains G the kernel (13) 
can be simplified. In particular, if the domain G has 
two axes of symmetry, and its upper and lower limits 
are described, respectively, by the functions y = +r\(x) 
at x e [-1,1] , (13) is real and it takes the form 



2n 



sin(c 2 (4-^ 2 )ri(x)) 
(s' 2 -s 2 ) 



dx . 



(14) 



In the case of a rectangular aperture (13) has the form 



^(*l j *2 ' *1 ' s 2 > c l ' c 2 ) — 



sinq^j' -Sj) sinc 2 (s 2 - s 2 ) 
n(s[ -Si) n(s' 2 -s 2 ) 



(15) 



Note, the Equations (10) and (12) are equivalent in 
this sense [3]: between solutions of these equations 
there exists bijection, that to every solution of (10) 
corresponds the solution of (12) and on the contrary. 
If /„ is the solution of (12), then corresponding 
solution U t of (10) is determined by the formula 

U t (x,y) = -^^F(Q)exp[i(wLgf t (Q)-(c l xs l +c 2 ys 2 ))]ds l ds 2 ;(16) 

if U„ is the solution of (10), then corresponding 
solution /„ of (12) is determined by (4). 

Using a more general expression (13) for K(Q,Q',c), 
we consider the operator 

Df = \\K{Q,Q\c)f{Q')dQ' (17) 

Q 

and corresponding to it quadratic form 

(Df,f) = \\\\K{Q,Q',c)f{Q')dQ'J{Q)dQ = 
n n 



= ^\\\ jjf(Q)ox P (ic(P,Q))dQ 
(2n) g ^ 



dxdy > . 



(18) 



where 



From here follows that K(Q, Q',c) is a positive kernel 
[7]. Accordingly the operator D is also positive on the 
nonnegative functions cone 7C of the space C(G) [8]. 
Based on this property the operator D retains 
invariant cone K , that is DK a 7C . 

Since a set of values of operator A is a set of 
continuous functions, belonging to the space L 2 (D.) , 
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and a set of continuous functions in the domain Q , is 
dense in the space Lj(Q) [7], we shall investigate the 

solutions of (12) in the complex space C(Q) . 

On the basis of decomplexification [9] we consider the 
complex space C(Q) as a direct sum 

C(Q) = C(Q) 9 C(Q) of two real spaces of continuous 
functions in the domain Q . The elements of this space 
have the form whose elements are given in the form: 
f = (u,v) T eC(D.), w=Re(/)eC(Q), v = Im(/) e C(Q) . 
Norms in these spaces we shall introduce as: 

I U llc(Q) = ™1 1 U ^ I ' I V llc(fl) = I I ' 

|/lc(fi) =max (l M lc(fi)'ll v lc(a))- ( 19 ) 

The Equation (12) in decomplexified space C(Q) we 
reduce to equivalent to it system of equations 

«(G) = B,(«,v) - JJf(6')^(6,8',c) M(g,) =48', 
a V M (2') + v 2 (e') 

v(G) = B 2 (".v) = |fF(«2')g(j2,8',c) v(g ' } rfg' . 
n V M (6') + v 2 (e') 

(20) 

Denote a closed convex set of continuous functions as 
S M a C(Q) supposing that 

= V © V • V = {" e V = ||"|| C(fl) * , 

S K ={veS M/ .\\v\\ cm <M}, 
M = max [ [ F(g') I K(Q,Q', c)\dQ . 

Consider one of the properties of the function 
exp(;arg /(£')) included in (12) at f(Q')^>0. It is 
obviously that 

(« 2 (G') + v 2 (Q')) 

is continuous function if u(Q') = Re f(Q') and 
v(G') = Im f(Q') are continuous functions, where 
|exp(iarg/G2'))| = l for any f(Q') . If u(Q')^>0 and 
v(Q') -> simultaneously, then /(£?') = is a complex 
zero, the argument of which is indefinite according to 
definition [10, p. 20]. On this basis we redefine 



exp(iarg/G2')) at u(Q')^Q and v(Q')^>0 as a 

function which module is equal to unit and indefinite 
argument. 

Theorem 1. The operator B = {B l ,B 2 ) T determined by (20) 
maps a closed convex set S M of the Banach space C(Q) in 
itself and it is completely continuous. 

To prove the theorem it is necessary to show that 
B : C(Q) C(Q) and using the Arzela Theorem [6] to 
prove that functions of a set S g = BS M are uniformly 

bounded and equapotentially continuous where the 
inequality BS M cz S M holds. 

As a corollary of Theorem 1 the satisfaction of the 
Schauder principle conditions [11] in accordance to 
which the operator B = (B l ,B 2 ) T has a fixed point 

/„ = (k, ,v f ) T , belonging to a set S M , follows. This 
point solves (20) and, respectively, (12). Substituting 
/* ={u*tV t ) T into (16), we obtain the solution of (10), 
which is a stationary point of functional (9). 

Concerning the synthesis problems of linear radiator 
for the case of one-dimensional domains Q , solutions 
of a system of equations similar to (20) were 
investigated, in particular, in [2, 3]. The obtained there 
results show that nonuniqueness and branching of 
solutions dependent on the size of one-dimensional 
physical parameter of the problem, are characteristic 
for equations of the type (12). The results [2, 3] can not 
be transferred directly to two-dimensional nonlinear 
integral equations (12), (20). Here, unlike the 
branching points [3], there exist branching lines of 
solutions, and the problem of finding the branching 
lines is a nonlinear two-parameter spectral problem. 

We shall formulate important properties of (12), 
necessary further, which are checked directly. 

1°. If the function f(Q) solves (12), then complex- 
conjugate function f(Q) is also a solution of (12). 

2°. If the function f(Q) solves (12), then e' T /(2) is 
also a solution of (12), where x is an arbitrary real 
constant. 

3°. For functions F(s l ,s 2 ) even on two arguments (or 
on one argument) linear operator B , that is in the 
right part of (12), is invariant concerning the type of 
parity function arg f(s l ,s 2 ) on two arguments (or 
concerning the argument on which the function 
F(s l ,s 2 ) is even). 
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Further for the uniqueness of desired solutions using 
the property 2° we shall select the parameter x to fulfill 
the equality 

arg/(0,0) = 0. (22) 

It is easy to check that for the case of a symmetric 
domain Q the function 

MQ,c) = \\F{Q')K{Q,Q',c)dQ' (23) 

is one of the solutions of (12). Further we shall call 
solution (23) as primary. Since, as shown above, the 
operator D defined by equality (17) is positive on the 
cone of nonnegative functions 7C e C(Q) and F a 7C , 
then / = DF is also nonnegative function in the 
domain Q . 

To find the branching lines and complex solutions of 
(12), which branch-off from the real solution f (Q,c) 
we consider the problem on finding such a set of 
values of parameters c* 0) = (cj -* , ) and all distinct 
from f (Q,c) solutions of (20), which satisfy the 
conditions 

max| M (e,c)-/(e,c (0) )|^0, max I v(Q, c) I -> (24) 

QeD. 1 1 QeQ 1 1 

as |c-c (0) |h>0. 

These conditions mean that it is necessary to find small 
continuous in Q solutions 

w(Q, c) = u(Q, c) - / (g, c (0) ) , (0(6, c) = v(fi, c) , 

which converge uniformly to zero at c— >c (0) . In 
addition, it is necessary also to take into account the 
direction of convergence of vector c to c (0) . 

Put 

q^c^+ii, c 2 =4 0) +v (25) 
and we shall find the desired solutions in the form 

«(e,c) = / (e,c (0) )+w(e,n,v), v(g,c)=co(g,^v). (26) 

Further we omit the dependence of functions 
w(Q, |i,v) and a>(2,|i,v) on parameters \i and v to 

simplify the notations. 

Indicate the properties of integrands in the system (20): 
they are continuous functions of arguments. 
Substituting (25) and (26) in (20), we expand 
integrands in the uniformly convergent power series 
with respect to the functional arguments on functional 



arguments w, a> and numerical parameters |i and v 
in the neighborhood of point (c (0) , /„(£), c (0) ),0) : 

F(Q')K(Q,Q',c) U(Q \ = 

= Z We>eV> m (e>"(<2>^ 9 , (27) 

m+n+p+q>0 

F(Q')K(Q, Q', c) V(Q '\ = 

V« 2 (e')+v 2 (e') 

= Z B mnpq {Q,Q\^y)w m {Q') b f{Q')v P v q . (28) 

m+n+ p+q>\ 

Here A mnpq {Q,Q' ,c (Q) ) , B mnp<! (Q,Q',c m ) are the 

coefficients of expansion continuously dependent on 
arguments. Substituting (25) - (27) in (20) and 

considering / (g',c <0) j as a solution of (20), we obtain 

a system of nonlinear equations with respect to small 
solutions w and o>: 

w(e)=fl 1 o(e,c (0) )n+ flol (e,c (0) )v+ 
+ x n /, v^jA ffl „ / , 9 (e,e',c (0) )w" i (e>"(eve' ,(29) 

m+n+p+q>2 Q 

.(Q)-\\F( Q) K(Q,Q^)-^^ d Q' = 

= Z n p v'j(B MM (ae',c (0) )/(fi>"(e')rfn^ , 

m+n+p+q>2 q 

(30) 

where 

fllo (e,c (0) )=jjA, 01 o(e,G',c (0) )rfG' 
flol (ac (0) )=jjA)ooi(e,G',c (0) )rfe'. 

n 

For further application of the methods of the 
branching theory of solutions of nonlinear equations 
[12] to system (29) and (30), it is necessary to find 
solutions distinct from the trivial of linear 
homogeneous integral equation. We obtain this 
equality equating to zero the left part of (30): 

9(2) = T( Cl ,c 2 )q> = jj f . F n iQ) K{Q, Q', c, ,c 2 ) cp(G') <*G' (31) 

at / (G',c)>0 . Note, that the operator 
T(c) : C(Q) -> C(Q) is completely continuous. 

Accordingly to [12] those values of parameters 
(cj 0) ,C2 0) ) e M 2 at which the linear homogeneous 
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equation (31) has solutions distinct from identical zero 
are points of possible branching of solutions of 
nonlinear equations (29) and (30). Eigenfunctions of (31) 
are used at construction of branching-off solutions of 
(29) and (30). 

Note, in a general case (31) is a nonlinear two- 
parameter spectral problem. 

Implicit Function Method for Finding the 
Connected Components of Spectrum of 
Nonlinear Two-Parameter Spectral 
Problems 

Here we shall consider a method for solving one class 
of nonlinear two-parameter spectral problems, which 
include, in particular, the problem on finding the 
solutions of linear homogeneous integral Equation (31). 
Since the proposed approach can be applied to finding 
not only the connected components of the spectrum of 
integral equations, we present the material on the 
general (operational) level. 

Let E and V be a complex Banach space and a vector 
parameter X = (X 1 ,X 2 ) belongs to the domain (open 
connected set) A = Aj x A 2 of the complex space 
C 2 = CxC . Here X i e A, c C , A,. = [X t e A,. : \X t \< r x } 
(i = 1, 2) , r x is some real constant. Consider the 
operator-function A (•,•) : A — >L (E,V) , where to each 
X = (X x ,X 2 )eA the operator A(X l7 X 2 ) eL(E,V) is 
assigned. Here L(E,V) is the space of linear bounded 
operators [9]. 

Let us consider the nonlinear two-parameter spectral 
problem of the form 

A(X l ,X 2 )x = 0, (32) 

in which it is necessary to find the eigenvalues 
X = (xf^ , X 2 Q) ) e A and corresponding to them 

eigenvectors x (0) e E ( x (0) * ) such 
thatA(^ \4 0) )x (0) =0. 

Without limiting the generality, in view on the form of 
(31), we shall present the operator-function A(X l ,X 2 ) 
as 

A(X l ,X 2 ) = T(X l ,X 2 )-I . (33) 

Here T(X l ,X 2 ) is a linear completely continuous 
operator acting in the Banach space E and analytically 
dependent on two-dimensional parameter (X i ,X 2 ) , I 



is a unique operator in E . In this case we put that 
A (•,•): A ->L (E,E), i.e. V =E, and A(X U X 2 ) sL(E,E) . 

Let the Banach spaces E and E n (n = l,2,...) and also 
the system P=(p„) of linear bounded operators 
p„:E -> E n such that 

Ikvik \\4e (neN) Vie£ (34) 

be given. Operators p n are called conjunctive 
operators [13]. From the principle of uniform 
boundedness for p n the inequality jp n || < const 

(n e N) follows. 

Let in every space E n the element x n be selected. 
Writing these elements in order to increase the 
numbers we shall form a sequence {x n } . 

Definition 1 [13]. The sequence {x„} neN , from x n e E n 
P converges to xeE, if \x n - p n x\ E — > (neN'); we 

p 

denote x n >x (neN' ). 

As noted above, in this definition the elements x and 
x n belong to various spaces. The properties of P - 
convergence are given, in particular, in [13]. At 
application of that or other approach to discretization 
of initial problem the operator-function A{X l ,X 2 ) is 

approximated, respectively, by the approximate 
operator- functions A n (X 1 ,X 2 ) , neN. As a result, for 

each X = (X l ,X 2 ) e A we obtain a sequence of 

operators A n eL(E n ,E n ) converging to operator 

AgL(E,E) at satisfaction of the corresponding 

conditions. We shall present the definition of 
convergence of the operators A n eL(E n ,V n ) to 
AeL(E,V) necessary later on. Let the Banach spaces 
E,V,E n ,V n , n = l,2,..., and linear systems P=(/?„) 
and Q=(q„) of linear (connected) operators 
p n :E^E n and q n : V -> V n such that 

\\p n x\\ -> ||jc|| and ||^ B y|| -> |y|| (neN ) Vx e E , y eV 

(35) 

be given. 

Definition 2 [13]. The sequence {Ai} neN of the operators 
A n :E n ^V n PQ- converges (or discretely converges) to 
the operator A : E -> V , if for each P - convergence 
sequence {x n } the following relation is valid: 
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->x(neN) => A n x n — — >Ax; 



(36) 



we denote A n — -^-> A («eN)orA B ^A (neN). 

Discretization of initial problem (32), the choice of the 
spaces E n and determination of operators p n :£—>£„, 
q n : V — > V n are realized by various ways. In particular, 
if the operator-function is described by (33) and E is 
the separable (infinite-dimensional) Hilbertian space, 
one of the approaches to discretization (32) consists in 
the following. Consider an arbitrary complete 

orthonormal in E system of functions • Each 

CO 

element xeE can be presented as a series x = ^c k x k . 

k=i 

Here c k =(x,x k ) is the Fourier coefficient of the 

element x . If T(X t ,X 2 ) is a linear continuous operator, 

acting in the separable Hilbertian space, then matrix 
representation has the form [14]: 



T M (X 1 ,X 2 )-[t jk (X 1 ,X 2 )). k=i . 



(37) 



result, the approximated operator T M (X) is 
described by the finite-dimensional matrix-function 



T M„ (^l'^2) ~(tjk (^l'^2)) j k=1 • 



(39) 



Applying other methods of discretization to problem 
(32), in particular, the quadrature (cubature) processes 
for the case of homogeneous integral equations and 
change of derivates by their difference analogues in 
differential equations, we obtain approximate 
problems to find approximately the eigenvalues and 
eigenfunctions in the form 



A n (X l ,X 2 )x n =0, neN. 



(40) 



Moreover, the problem of determination the 
eigenvalues is reduced to finding the roots of the 
determinant of n -th order that is the roots of equation 

a n (X i ,X 2 ) a l 2(X- l ,X 2 ) ■■■ a ln (X l ,X 2 )^ 
a 2l (X l ,X 2 ) a 22 {X^,X 2 ) ••• a 2«(^i>^2) 



¥ n (^ 2 ) = det 



y a nl (X u X 2 ) a n2 {X x ,X 2 ) ... a m {X x ,X 2 ) 



= 



Here tj k (X l ,X 2 ) = \T(X l ,X 2 )x k ,Xj\ . This sequence of 
Fourier coefficients of the element y = T(X x ,X 2 )x is 

derived from the sequence of Fourier coefficients of the 
element x by means of transformation by the 
matrix T M (X V X 2 ) ■ 

Using the matrix representation of operator T M (X l ,X 2 ) 
the spectral problem (32) is formulated as 



A M (X 1 ,X 2 )x = (T M (X 1 ,X 2 )-I M )x = / 



(38) 



where I M is a unit matrix in the space of sequences 
l 2 .Thus, the operators T(X) and T M (X) are equivalent 

in the sense that they assign one and the same element 
y e E to the same element x e E . But we obtain the 
Fourier coefficients of element y = T(X)x as a result of 
action of the operator T M (X) on the element x . 

Obviously, that the spectrums of these operators 
coincide, that is the spectral problems (32) and (38) are 
equivalent. 

In this case we shall put that the finite dimensional 
spaces E n are generated by the bases {-^i}^ (neN) 
and to each element x e E the operators p n : E -> E n 

n 

assign the element x = ^c k x k , where c k =(x,x k ) . As a 

k=l 



(neN). (41) 

Note, if the coefficients a i j(X l ,X 2 ) are continuously 

differentiable functions of arguments, the partial 
derivatives of the function ^¥ n ( y X l ,X 2 ) are found by 

the rules of determinant derivation: 



8k i 

a 2 i(X|,A, 2 ) 



=2> 



BX, 
dk, 



dk, 



j = 1,2. 



fl lJ+l(^l>^2) 



n(\'^-2 



(42) 



Consider the auxiliary nonlinear one-parameter 
spectral problem, necessary later on, as a special case 
of (32). Assume that variable X 2 in the operator- 
function A(X 1 ,X 2 ) is expressed by some one- valued 
differentiated function X 2 =z(X 1 ) mapping the 
subdomain A ljP c: A, into some subdomain A 2 p cz A 2 . 
In the simplest case we put X 2 = $X t ( p is some real 
parameter). Introduce the operator- function 
Ap(A, x ) sA^.z^,)) at X,eA ip (reduction of the 
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operator-function A(X,,X 2 )) with which is connected 
nonlinear one-dimensional spectral problem 

A ? (k 1 )x = 0. (43) 

Here we assign to each value % = (X 1 ,z(X 1 )) e A the 
operator i4 p (^ 1 ,z(A, 1 )) eL(E,E) . 

Analogously to (40), we consider an approximate 
sequence of discrete problem of (43) 

Ap, n (Vz(^)k=0, neN. (44) 

Denote the spectrum of operator-function A p (X[) as 
s(Ap) . Assume that s(A^) * A ip . For the spectrum 
s(A ) of the problem (34) the following theorem is valid. 

Theorem 2. Let the following conditions be satisfied: 

1) operator-function A( v ) : A ->el_ (E,E) is holomorphic, 
and s(A ) * A ; 

2) operator-functions A n (•,•): A ->eL(£ n ,£ n ) are /zo- 
lomorphic and for any closed bounded set A c: A ilze 
following inequality max\\A n (X l ,X 2 )\\< 



< c( A ) = const (n e N) is calid; 



3) operators A(X U X 2 ) gL(E,E) , A n (X,,X 2 )e eL(E n ,E n ) 

( n e N ) are ffte Fredholm operators with zero index for 
any X = (X X ,X 2 ) e A ; 

4) spectrum s(A p )^A ip and a sequence of functions 
x ¥ n (X l ,X 2 ) are differentiable in the domain A; 

5) A n (X)^A(X) is stable for any X e r(A) = AU(A) . 
Then the following statements are true: 

1) every point of spectrum X{ 0) e.y(A p ) is isolated, it is 
eigenvalue of the operator i4p(X, 1 )=A(X, 1 ,z(A, 1 )) , the 

finite-dimensional eigensubspace Af|A(xj 0) )j and the 

finite-dimensional root subspace correspond to it; 

2) for each X[ 0> e s(A p ) there exists a sequence {X^J from 

X[ ^ e s(A M ) (n > n ) , swc/i i/iai X^ -> X,? ; 

3) eac/i poini X <0) = |x[ 0) ,z(x{ 0> jj e A is a spectrum point 
of the operator-function A (X, , X 2 ) ; 

4) if in some smaii s - neighborhood of the point 
X <0) =(x[ 0) ,z(x i (0) ))e A for all n, larger than some 
number N Q , the sequence of partial derivatives 



^^(^i,«>z(^i,n))j iS nonzero, then in an arbitrarily 

small e, -neighborhood of the point (x{°\z(x< 0) ))eA 

there exists a continuously differentiable function 
^2,n, =( Pw. (^i) which is a solution of (41), where 

Xi ^ = cp^ |x{ ^, j and at the point 

=( x( "k^N,( x ik)) u arbitrar y smal1 

deviates from the point spectrum of auxiliary one- 
parameter problem (43) |xj 

some 



(0) 

N, 



,(0) 



That is, in 

bicylindrical domain 

A ={(^1^2) e A : |^l-^l 0) | <e l' |^2-4 0) | <e 2} 

there exists a connected component of spectrum of the 
operator-function A W> (X 1 ,X 2 ) ( e l , e 2 are small real 

cons Ian Is). 

Proof. The proof of Theorem is based on Theorems 1, 
2 [13, pp. 68, 69] and on existence of implicit functions 
(see, for example, [15]). At first we show that the 
conditions of Theorem 1 [13, p. 68] concerning the 
existence of discrete spectrum of operator-function 
A p (Xj), follow from the conditions of formulated 

Theorem. Under the conditions of Theorem the 
operator A(X P X 2 ) is Fredholm operator with zero 

index for each X, = (X P X 2 ) e A , and the operator- 
function A (•,•) : A ->el_ (£,£) is holomorphic. From 
here follows that for each X, e A l p the operator 
A p (Xj), as reduction of the operator A (Xj,X 2 ) , is also 

Fredholm operator with zero index, and the operator- 
function A p (X|) : A] p ->L(E,E) is holomorphic. So, 

for the operator-function A p (X l ) the conditions of 
Theorem 1 [13] are satisfied. From Theorem 1 follows: 
each point X, (0) e s( A p ) is isolated, it is the eigenvalue 
of the operator A p (X,) , the finite-dimensional 

eigensubspace and finite dimensional root subspace 
corresponds to it. Thus, each point 

X (0) = (x< 0) ,z(x< 0) ))eA is the spectrum point of the 
operator- function A(Xj,X 2 ) . 

In addition, the conditions of Theorem 2 [13] are 
satisfied for the operator-function 

A p (X,): A ip — >L(E,E) . From this theorem follows: at 

n larger than some n e N for each X[ 0) e s{A p ) there 

exists a sequence {X ln | such that X [ n e s(A^ n ) and 

X U ^X[ 0) as «h>oo. Thus, each point (x< 0) ,X 2 0) ) = 
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= \X\ ,z[X\ 'jj is the eigenvalue of operator-function 
Ap(^j): A ip -*L(E,E) and, respectively, the 
eigenvalue of operator-function A (•,•) : A ->el_(£,£) . 
Since X[ ^ es(A pn ) is the root of (41), then from here 
follows that for n e N 

V n , z (*£ )) - <P„ . z ( *f )) = as « -> co . (45) 

From the criterion of convergence of sequence 
^°'ei(A M ) to X,{ 0) es(A p ) follows that for an 
arbitrarily small number e, > there exists such 
number N„>n that j^j ^ -A,[ 0) |<e, and 

^.(C^(C)) = o- 

Let X 1 , X 2 be independent variables in the domain A , 
and (x{ 0) ,4 0) ) = (A.{ 0) ,z(A.{ 0) ))eA be a spectrum point 
of the operator A(A,j,A, 2 ) belonging to s(A p ) cz s(A) . 

Under the conditions of theorem the functions 
x ¥ n (X l ,X 2 ) are differentiable in the neighborhood of 

3X1/ 

the point (x[°\X 2 0) ) and ^fe, z (^)Wo. In 

addition the point (x{ ^ , z(^{ ^ )) belongs to e, - 

vicinity of the point (x{ 0) ,X 2 0> ) . According to the 
Theorem about implicit function [15] in some 
neighborhood of point fa 0) ,\f>) there exists the 
continuous differentiable function X 2 =(p N ^(X 1 ) , 
solving the equation (41), where <p n {xf > ^ ) = z(^°^ )■ 

From here follows the existence of connected 
component of spectrum of the operator-function 
A(; •) : A — >L (E,V ) in some bicylindrical domain 

A ={(X 1 ,X 2 )eA :\X 1 -^ 0) \<e 1 ,\X 2 -Xf ) \<s 2 \. 

At the point (C.'^.) = (C'^,(^.)) this 
component deviates from the point of spectrum 
X (0) =(^,z(i! 0) ))eA of operator-function A(^,A, 2 ) 
arbitrarily. Theorem is proved. 

Consider Equation (41), putting that the coefficients 
a y (A,j,A, 2 ) of the matrix-function A B (A,j,X, 2 ) are 

differentiable continuous functions, i. e. partial 
derivatives d"¥ n (X 1 ,X 2 )/dX 1 and ST n (X l ,X 2 )/dX 2 
determined by (42), exist and they are continuous. If at 



the points {X^X^} e A i xA 2 the derivative dWjdX 2 
is nonzero, then solving the Cauchy problem [16] 

dX 2 _ d^JX^/DX, 

dx x w n {X v x 2 )idx 2 ' 

X 2 (X\ V> ) = X[ V} , (47) 

in the vicinity of each point X\ v) e A, we find v -th 
connected component of the spectrum of the operator 
function A n (X { ,X 2 ) . 

Numerical Finding of the Branching Lines of 
Solutions 

Return to finding the solutions of (31), where c, , c 2 
are spectral parameters. Let (c, , c 2 ) e A c , 
K = A q xA C2 , where 

K, ={ c « eA c, :0<c,.<r r }. 

We shall present equivalent to the previous definition 
of eigenvalues and eigenvectors for (31), which later 
on will be used for finding the numerical solutions of 
this equation. 

Definition 3. We shall call a set of values of parameters 
(q (0) , c 2 0) ) e A c at which n = 1 

|J(P(G) = If ffn^ ^ * (2, 2', q , c 2 ) cp(G W , (48) 

as eigenvalues of (31), and distinct from identical zero 
solutions that correspond to the values of parameters 

c} ' , c 2 0) J e A c are called eigenfunctions. 

By direct check we set that for arbitrary values of 
parameters (q,c 2 )eA c one of the eigenfunctions is 

q>o(G,c) = \\F{Q')K(Q,Q\c)dQ' = f (Q,c) . (49) 
n 

Write, necessary later on, conjugate with (31) equation 
V(0 = (c) V = -f^- f J K(Q, G'.c)v(G') 42' • (50) 

At arbitrary (q,c 2 ) e A e one of the eigenfunctions (47) 
is 

V Q (Q) = F(Q) (51) 

that unlike (pj(g,c) does not depend on q, c 2 . 

The existence distinct from identical zero solutions of 
(31) for arbitrary (q,c 2 ) e A e indicates the existence of 



( 
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connected components of the spectrum, that coincide 
with the domain A c . So, the condition of Theorem 2 is 

not satisfied: s(A ) * A c . To fulfill this condition we 
eliminate the eigenfunctions (46) from the kernel of 
(31), namely, we consider the equation 



9(6, c) = f(c)q> = jj K(Q, Q\c)m') dQ' , 



(52) 



where 



K(Q,Q',c) = 4^K{Q,Q\c)- ^f Q ) C) . (53) 



/o(G',c) 



<Po 



From the Schmidt Lemma [12, p. 132] follows that 
u = l for any values (q,c 2 ) Equation (53) not be an 
eigenvalue in the sense of Definition 3, that is q> Q (Q,c) 
will not be an eigenfunction of this equation. The 
connected component, which coincides with the 
domain A c and corresponds to function (p (g,c) , is 
excluded from the spectrum of operator. 

On the basis of (53) we are convinced that the 
following inequality is valid 



<+oo (54) 



qc 2 


c,c 2 


(2k) 2 


_(2,) 2 



for kernel of operator f (c), where 
F(Q) 



a, 



max 

ceA, 



/o(G,c) 



a. 



JJi 



F(Q) 



/o(G,c) 



-dQ, 



ff |/o(6',c)| 
a / •• ii <*G • 



J dl|/o(2',c)| 



From inequality (54) follows that the operator T(c) is 
Fredholm with zero index [7]. Moreover, it is 
completely continuous operator in the space (Q) . 

Functions included in the kernel of (52) admit an 
analytic continuation in the complex domain A c , if q , 
c 2 is considered as complex parameters. Holomorphy 
of operator-function A (q , c 2 ) = f (q , c 2 ) - / follows 
from the differentiation of the integral operator 
r(q,q>) on parameters q , c 2 ■ In particular, the 
existence of partial derivatives 3A(q,q,)/3q , z'=l,2, 

at any point (c{ 0) , c 2 0) j e A c , follows from the 

continuity of the kernel K(Q,Q',c) on the set of 
variables in the domain QxQxA c and the existence 
and continuity of the partial derivatives 



SK(Q,Q', c^c^/dc^ i = 1, 2 , in Q x Q x A c . It is easily to 
be convinced in this. 

Putting in (52) c 2 = (3q ( P is a real constant) we 
consider a nonlinear one-dimensional spectral 
problem 

9(G) = T(q)9 - Jj£(G,G',<i>p(GW ■ (55) 

Here k(Q,Q',c l ) = K(Q,Q',c l ,$c l ) . Since f(q) is 
narrowing of operator T(c) , from here follows that 

T(c) is Fredholm operator with zero index at 
arbitrary c, sA, , and the operator-function 

A (■)= (f( •)-/): Ai (E,E) is holomorphic. If 
i(A ) ■£ Aj , from holomorphy of operator-function and 

Fredholmic of the kernel K(Q,Q', Cl ) = K(Q,Q', Cl ,$ Cl ) 
follows satisfaction of conditions of Theorem 2 
according to which each point Cj (0) e s(A ) is isolated 
and is eigenvalue of (55). Respectively, the points 
(cfV^) = {4 0) >& c i 0) ) are eigenvalues of (52). For the 

approximate finding the connected components of the 
spectrum of (52) in the vicinities of the points 
(cf^c^j we solve the Cauchy problem (46), (47), 

where as the initial conditions the found solutions of 
auxiliary problem |cj (0) , c 2 0) j = |cj <0) , Pc} ' j are used. 

Pass to construction of algorithms for numerical 
finding the eigenvalues of (31) and (52), using the 
discretization methods [9, 13]. 

Consider any convergent cubature process [13]: 
\\x(Q)dQ = f j a jn x{Q jn ) + ^ n ( X ) (neN) (56) 

with the coefficients a jn e K and nodes Q jn e Q , 
j = 1 -s- n . Reject the remainder term § n (x) in (56) and 
replace integral in (55) by it. Giving the variable Q 
the values Q = Q in ( i = \+n ), we have the 
homogeneous system of linear algebraic equations 
concerning u ln ,...,u nn \u ln ,...,u nn : 

n 

U in = Tm„ ( C l )"„ = Z a jn K (Qin . 2j« . C l ) « jh ( ' = L « )/ (57) 
7=1 

where u in =u(QJ. 

Consider the question of solvability of system (57) 
and convergence of its solutions to solutions of (55). 
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For this purpose, we consider (55) and the system of 
equations (57) as operator equations in approaching 
spaces. The integral equation (55) we shall consider in 

the Banach space E = C(Q) , where f (q) eL(E,E) . 

We consider the system of linear equations (57) as 
opera- 



tor equation 

u n=f M M)u n (neN) 



(58) 



in the Banach space £„=C(Q„) , consisting of 
vectors — grid functions u n = (u ln ,u 2n ,...,u nn ) . Here 



IMk = F»llc(n„)=™^l M » 



and operator f M el_ (E n ,E n ) . The connected operators 
p n eL{E,E n ) are defined by the formula 

p n u = (u(s ln ),u(s 2n ),...,u(s nn )) for all u e E = C(Q) . 

Thus, the problem of the eigenvalues 
Ap(q)u = |f (q)-/ju = is approximated by the 

problem of the form A p ; „(c 1 )u„ = (f Mn (q) -/„ ju„ = , 
where A p „Q: A, -+L(E n ,E n ) . 

According to [13, p. 99] for operators defined by (55) 
and (57) the following result is valid: if the 
kernel K : Q x Q — > R is continuous and the cubature 

process (56) converges, then T M — > T is the compact 1 . 

Since in these conditions the operators 



1 Comment 1. Note, for nonlinear one-parameter with respect to 
IeAcC problem of eigenvalues 

(i) x(t) = ]K(t,s,X)x(s)ds (A(X):=x-T(X)x = 0) 

D 

and corresponding discrete problem 

n 

(») *,„( f ) = Z fl > K (^'V^) x ;»' »' = 1.2,...,ra, 

7=1 

(A n (X)x n :=x n -T(X)x„=0) 
the condition 1-4 of Theorem 2 will be satisfied if the 
corresponding quadrature process is convergence and the kernel of 
the integral operator K (f, s,X) satisfies the following conditions: 

K (t,s,X) is Holomorphic relative to XeAcC; 

K(t,s,X) and dK(t,s,X)/dX are continuous relative to 

(t,s)eDxD; 

K (t,s,X + AX) -K (t,s,X) dK (t,s,X) 



max 

reD 



AX 



ex 



ds -> at AX -> . 



f( Cl )el_(C(Q),C(Q)) and f Mt (X) eL(c(Q„),c(Q„)) 

are completely continuous, then the conditions of 
Theorem 2 are satisfied and on this basis we have: 



1) for every c[ 0) e s(A) there exists a sequence 
{c^}ss(A n ) (n>n Q ) and c<°> -> q (0) ; 

2) if q ( n 0) ^ q (0) e Aj , {c™ ) eJ (A„) (neN'czN), 
then c, (0) £s(A B ). 



Thus, solving the eigenvalue problem (57) in a finite- 
dimensional spaces E n = C(Q„), we find approximate 
eigenvalues that converge to the exact solutions of (55) 
at n — > oo . Finding the eigenvalues of (57) is reduced, 
in particular, to finding the roots of the equation 



*„(c 1 ) = det(f Mn (c 1 )-/„) = 0. 



(59) 



In addition, it is necessary to introduce such condition: there exists 
such X e A which is not eigenvalue of integral equation (i) [13]. 



Return to two-dimensional spectral problem (52). 
Applying (56) to (52), we have a system of linear 
algebraic equations 

n 

U in = T M„ ( C l ' C 2 ) U = E °jn ^ (Qin > Qjn > C l > C 2 )" jn (60) 
7=1 

similar to (57). 

Consider finding the solutions of the equation 

W n ( Cl ,c 2 )=det(T Mn ( Cl ,c 2 )-I n ) = (61) 

as a problem on finding the implicitly given function 
c 2 =y(c 1 ), reducing it to Cauchy problem (46) and 
(47). We use solutions of (59) as the initial conditions 
(47) since to each isolated root of this equation cff 

corresponds eigenvalue 
cT=(cS ) ,c™) = ( e g>,z(cg>)) e A e of problem (60). 

Thus, the initial conditions (50) in the Cauchy 
problem we determine as cfj = fic^ . If 

BW i \ 

" \ c U' c 2i K ' we solve the problem (46), (47) in 

dc 2 ^ ' ' ' 

each vicinity of points (c^.c^) = (c 1 ( ° ) ,Pc 1 ( ° ) ) e A c . 
Then we find continuously differentiable function 
y,- (q) that satisfies the condition y( c i?) = 47 • 

In the case when c[ 0) =(q ( ° ) ,Pq ( ° ) )e A c is a real 

eigenvalue of the problem (59), y, (q ) are a real 
differentiable functions that describe some smooth 
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curves in the vicinity of points c^ 0) = (q ( °\Pq (< f) e A c . 

That is, in this case (60) has a line spectrum. 

Solving the problems (59), (60), we find a set of 
parameter values (q,q,)eA c at which from the real 

initial solution / (g,c) the branching-off of complex 

solutions of (20) is possible, in which functional (9) 
takes smaller values than the solution f (Q,c) ■ 

Show the basic relations for reduction of two- 
dimensional integral equation of (31) to the 
corresponding system of linear algebraic equations of 
the type (60) using the Gauss quadrature formulas that 
are necessary at calculations by means of the libraries 
of standard programs. 

Put that Gauss formula has m nodes. We introduce 
into consideration the vector of dimension n = m 2 

U = \u xx ,U x2 ,...,U Xm ,U 2x ,U 22 ,...,U 2m , ...,M ml ,M m2 ,...,M mm j = 

= {u(s 1 ,s 1 ),u(s x ,s 2 ),...,u(s 1 ,s m ),u(s 2 ,s l ),u(s 2 ,s 2 ),...,u(s 2 ,s m ), ... 
•••> "Om . h)Ms m , s 2 ),..., u{s m ,s m )} = [u{ Si , Sj)}" =i = \u Vj }" , =i 



(62) 



Reduce four-dimensional matrix to the two- 
dimensional one with coefficients of the form 



( 



Yj°q 77 - P ' — t K ( S k> S l> S p' S q> C l> C 2) 
9=1 J0\ S p' S q' C \' C 2) 

p,q = \^m, k,l=\^m 



where b p ,b q are the Gauss nodes. 

Note, that solving the Cauchy problem (46), (47) by 
Runge-Kutta or Adams methods it is necessary to 
calculate the derivatives of function x ¥ n (c l ,c 2 ) at 

points of given grid. Such calculations can be 
performed by two ways. The first of them consists in 
direct differentiation of the determinant by (42). The 
second way is to find approximate derivatives using, in 
particular, the difference formula [17]: 



2h 



Taking into account the form of right part of (46) to 
find the derivative at the point (q,q>) we obtain the 
following expression 



dq h 

dV n {cj,4)~ V n {cj,cf l )-V n {cl,cf l ) 
dc 2 h 

_ ^ n {fi\c{)-^ n (c^,ci) 
"^(cl.cf 1 )-^^,^- 1 )' 



63) 



Consider the application of this method to numerical 
finding the connected components of the spectrum 
(spectral lines) (31) for concrete given 
functions F(s l ,s 2 ) ■ 

To find the spectral lines of (52) by solving the 
Cauchy problem (46), (47) it is necessary to find the 
initial conditions (47), i. e. to find the roots of (59). For 
finding the approximate solutions of this equation 
using the definition 3, on the basis of (55), we find 
approximately the set of parameter values q at which 
|i = 1 is an eigenvalue of equation in the sence of 
linear eigenvalue problem. Since in the synthesis 
problems of flat aperture the parameters q, c 2 are 
real and positive, we define the domain A e as: 

A c =|(q,c 2 ): 0<c l <d l , 0<c 2 <d 2 }. 

We shall find the approximation to the eigenvalues of 
(59) on the segments of beams c 2 = Pq belonging to 

A c as a discrete set of parameter values (q 7 -,Pq A at 

which a unit will be an eigenvalue of the matrix (57). 
For this purpose we select any mesh c Y - on the 

interval < q < d x . Then, in particular, using 
Danilevskiy method [18] we find the eigenvalues 
\Xj [c X j j , those values of parameter c 1; - at which 

\ij (q j ) « 1 will be approximations to the roots of (59). 

The next step consists in finding the corrected roots of 
(59), considering it as a transcendental equation with 
respect to q . In addition we use the solutions 
obtained in the preceding paragraph as initial 
approximation. We use the known methods, in 
particular, the Muller method [19] to find the roots of 
(59). As a result, we obtain the corrected values of the 

eigenvalues of (31) (q (0 j,4^) = (q (0 j,pq (0 J ?) on the 
beams c 2 = Pq that belong to A c . 

Consider a numerical example of finding the set of 
eigenvalues \ij (q <0 j) by Danilevskiy method. In Fig. 2 
the curves corresponding to variuos eigenvalues 
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on the beam c 2 =0.8cj for F(s l ,s 2 ) = l are 
shown. Vertical lines illustrate the parameter values 
c{°j for \ij (c[j j « 1 . The approximation to the 

eigenvalues of the problem (60) is given by bold points 
in Fig. 3. Approximate and corrected values of the 
roots of (59) at p = 0.8 are presented in Table 1. 

6:3 
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TABLE 1 



r (j) 


Approxim 
ate Value 


Corrected 
Value 


rP 


3.14 


3.1415927 


rP 


3.65 


3.6439078 


(3) 


3.92 


3.9269908 




4.88 


4.8679838 


r< 5 > 


5.35 


5.347 




6.28 


6.2831855 



Using the corrected values of the roots of (59) as initial 
data in the Cauchy problem (46), (47) for (60), we find 
the possible branching line of solutions of (20) which 
are shown in Fig. 4 for the function F(s x ,s 2 ) = 1 . 

Note, in the case of rectangular domain Q the 
variables can be separated [3] in (31) and 
F(s l ,s 2 ) = F^s^- F 2 (s 2 ) . In addition taking into account 
the form of kernel (15) and (23), we obtain two 
independent equations 




(64) 

They describe the possible branching line of solutions 
of (12), that form a rectangular grid in the domain A c . 

For function F(s 1 ,s 2 ) = 1 it is a grid of lines, which is 
shown in Fig. 4. In the case of nonlinear entering of 
spectral parameters into the kernel of integral 
operator (31), obtained by separation of variables of 
(64), do not describe all solutions of the initial 
equation (31). This follows, in particular, from 
analysis of Fig. 3 and Fig. 4. On the last of them we 
see that besides the direct spectral lines there exist 
also curves of the spectrum obtained first by solving 
(31) using the Cauchy problem. 

For the case when the variables in the function 



F(s l ,s 2 )= 



0, 



(65) 



FIG. 3 



are not separated, the spectral lines of (31) are shown 
in Fig. 5,b. In Fig. 5, a the solutions of auxiliary one- 
parametrical spectral problem (43) on different beams 
c 2 = Pq are given. These solutions are the roots of (59) 
and they are used for solving the Cauchy problem 
(46), (47) to find the spectral lines presented in Fig. 5,b. 

About branching of solutions of synthesis 
problem of flat rectangular aperture 
Partial cases 

Omitting the cumber some intermediate calculations, 
we shall present the main results obtained in [20, 21] 
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c 2 =Pc, 



FIG. 5, A 



C 2 =Cl 




FIG. 5, B 

when we studied branching of solutions of nonlinear 
two-dimensional integral equations (20) for the case of 
a flat rectangular aperture. We put also that function 
F(s l ,s 2 ) is even on both arguments, the domain Q is a 
rectangle and the multiplicity of the eigenvalue 
fcj" , cjP j e A c is two on the corresponding spectral 
lines, where 

A c ={(c,, c 2 )e A c : 0<c 1 <d l ,0<c 2 <d 2 }- 

The investigations of branching of solutions were 
carried out on the beams c 2 = Pq at < P < 1 . 

As noted previously, the function 

is one of the eigenf unctions of (31) for arbitrary values 
of the parameters (q, c 2 ) e A c . 

Except this function of eigenvalues 
c (/) =(c 1 (/) ,c 2 /, ) = (c 1 (/) ,pc 1 (/) ) of (31) belonging to the 

found spectral lines (see Fig. 4 and Fig. 5) there exist 
eigenfunctions distinct from cp . For the case when the 
multiplicity of an eigenvalue is two, we denote the 
second eigenfunction as <p 1 \s 1 ,s 2 ,c®,$c®) . In 



particular, for the first branching point , O.Scf 1 ' j 

(see Fig. 4 and Table 1) at F(s l ,s 2 ) = l the second 
eigenfunction has the form 

<p, (J, , s 2 , q (1) , pcf ) = s 2 ■ / (J, , s 2 , q (1) , pq (1) ) . (66) 

Giving analogously to (25) to parameter q (/) the small 
disturbance, i. e. putting 

Cj = C[ (/) c 2 = Pc[ (,) + P|i 

and taking for (20) similar transformation, we obtain a 
system of nonlinear integral equations of the 
Liapunov-Schmidt type concerning small solutions w , 
(0, similar to (29) and (30): 



w(Q)=a 01 (e,c ( ' ) V + 

X ^ }} A mnp 02, e',c (,) )w m osv (cw , 



(67) 



m+n+p>2 q 



«(G)-jj^(e')^(e,e',c (/) ) 



<B(<2') 



/o(G',c 



rf!2' = 



(68) 



= I n"j|B„ ffl;) (e,e'c ( 'V"(ev(eve', 

m+n+p^2 q 

where 

a 01 (e,c (0, ) = jjAo 01 (G,e',c (,) ),/G'. 



The left part contains two-dimensional linear integral 
equation of the type (31). Taking on the system 
transformations on the basis of general branching 
theory of solutions [12], we obtain that two complex 
conjugate solutions branch-off from initial (real) 
solution (23) at the point (c^tfcpA . In the first 

approximation they have the form 

f$ (s 1 ,s 2 ,c 1 ,Pc 1 ) = f (s 1 ,s 2 , c[ l) , p Cl (,) ) + 

+ (a ( ff , , s 2 , c\ l) , p C[ (/) ) + o$ (s v s 2 , c\ l) , p Cl ( '> ) hi ) jx ± 



±i 



(^,%,c{'>,pc«) 



V /2 +0( M 3/2 ). 



(69) 



Indicate some characteristics of branching-off 
solutions of (12) for the case of symmetric aperture G 
and even on both arguments (or on one argument) of 
the function F{s x ,s 2 ). In the first approximation the 
properties of solutions follow from (66). They are 
defined by the feature of functions 

/o^.ef.Pcf ), a(s l ,s 2 4\^ ) ), 
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a 02o( 5 'i' 5 '2> c i ,P c i J and the property of the 
eigenfunction cpj ^s l ,s 2 ,c{ l \$cl l) j which corresponds to 
the / -th point of branching of solution ^c[ l \^c^ j . The 

last two solutions have the same type of parity as the 
function f Q ^s l ,s 2 ,c ( l l) ,fic{ l) j . It is confirmed also by the 

general properties 1° - 3° of the solutions of (12) given 
previously in II. 

Numerical Algorithms of Finding the Optimal 
Solutions of Base Equations of Synthesis 

Give one of the iterative processes for numerical 
finding the solutions of (20), which is based on the 
method of successive approximations [3]: 

« B+1 (0 = B^vJ^jjK&Q'^FiQ') U " (Q '\ dQ' 
" K(Q') + v 2 n (Q') 



v„ +1 (0 = B 2 (u n ,v n ) = j\K(Q,Q',c)F(Q') V " (Q '\ dQ' 

(n = 0,l,...). (70) 

We substitute in (10) functions arg/ n (s p i 2 ) = 
= arctg(v n (s 1 ,s 2 )/« n (i 1 ,5 2 )) obtained through successive 
approximations (70). Thus we obtain a sequence of 
functions 

Un(Xt y)=Vl- ff F(S U S2)e {^fM-(my S2 )] (71) 

(2«) n 

which is denoted as {U n } . For a sequence {£/„} the 
theorem 4.2.1 with [3] is valid. It indicates that the 
sequence {U n } is a relaxational for (7), and numerical 
sequence {a(t/ n )} is convergent. 

In the case of even on both (or on one) arguments 
function F(s v s 2 ) and symmetric domains G and Q it 
is expedient to use the property of invariance of 
integral operators B l (u, v) and B 2 (u, v) in (20) relatively 
the type of parity functions u(s l ,s 2 ) , v(s v s 2 ) for the 
iteration process (41). Functions u , v having some 
type of parity on the corresponding argument belong 
to invariant sets O jjr O k( of space C(Q). Here indices 

i,j,k,£ take the values or 1. In particular, if 
u(s 1 ,s 2 ) e <D 01 then u(-s 1 ,s 2 ) = u(s 1 ,s 2 ) , and 
u(s 1 ,-s 2 ) = -u(s l ,s 2 ) . By direct check we make sure 
that there are such inclusions: 



B,(© s U©„)c:© #/ B 2 (3» s UOjc$ t( , 

From these relations follows the possibility of 
existence of fixed points of the operator B belonging 
to the corresponding invariant set, that is to the 
solutions of (20) and, respectively, to (10) and (12). 

Analysis of Numerical Results 

First we consider the analysis of branching of 
solutions of (12) with the kernel (15) for the given DP 
F(s l ,s 2 ) = l on the beam c 2 =0.8c 1 . The lines of 
possible branching of solutions were presented in Fig. 
4. The first few points of branching of solutions on the 
ray c 2 =0.8q are given in Table 1. From Fig. 4 and 

Table 1 follows that the point , Pq 0) ) = (71,0.871) is 

the first branching point of the initial solution 
/o( 5 i' 5 2)- Th e eigenfunction of linear homogeneous 
integral equation (31) 

<p t (s v s 2 , 71, 0.871) = s 2 ■ f Q (s 1 ,s 2 ,n,0.8n) different from 
(p and corresponding to that point is even on the 
argument Sj and odd on the argument s 2 . 
Accordingly branched-off at point 

Ci°\pci 0) J = (71,0.871) two complex conjugate 
solutions have the following properties: 
Re//^ (s l ,s 2 ,cl°\$cl° ) j is even on both arguments, 

and (j 1 ,J2» c i° ) »P c i° ) ) i s °dd function on 

argument s 2 . Respectively the functions 

arg/j^ (5 , 1 ,5 , 2 ,c}°\pcj ^ J are odd concerning s 2 . In 

particular, it is confirmed by the results of numerical 
experiment. 

The distribution function of oscillating velocity 
U(x, y) on the surface of the antenna which is 
determined by formula (71) depends on the properties 
of functions F(s u s 2 ) and arg f(s 1 ,s 2 ) • In the case of 

symmetric domains G and Q , even F(s x ,s 2 ) and 
function argf(s 1 ,s 2 ) odd on one of arguments, 
U (x, y) is a real and nonsymmetric function in respect 
to one of the coordinate axes of aperture G . 

Fig. 6 presents the value of the functional a F , which it 
takes on the primary (curve 1) and branching-off 
solutions (curves 2, 3, 4). Curve 2 corresponds to 
branching-off solutions with phase DP odd on s 2 . 
Optimal distribution function of oscillating velocity in 



( 
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aperture corresponding to this type of solutions is real 
and asymmetric relatively the axis OX. However, 
formed by it synthesized amplitude DP is symmetric. 

logc3 F 

y / 



~7 



zVi. 



FIG. 6 THE VALUES OF FUNCTIONAL ON VARIOUS TYPES OF 
SOLUTIONS 

Curve 3 corresponds to branching-off solutions at the 
point U 2) ,^ 2) ) = (3.6439,2.91512) with the following 

properties: 

\ w % fif Hi ,s 2 ) = -aig f$ {s u s z ), 
U^fif(si-s 2 ) = -^gf l f(s l ,s 2 ). 



(72) 



Eigenfunctions of (31) corresponding to this branching 
point have the form: 

9o(*i>* 2 ) = /o(*i>*2>3-6439,2.91512), 

q>! 0,,s 2 ) = s l s 2 f (s 1 ,s 2 , 3.6439,2.91512) . 

Curve 4 corresponds to the branching-off complex 
conjugate solutions at the point 

fcp.O.Sc} 5 ^ (5.347,4.2776) with the following 
properties: 

Jarg fif = arg fif ( s i> s 2)> 

arg fif {si~s 2 ) = arg fif ( s i> s i)> 



(73) 



that is, the optimal phase DP is even on both 
arguments. 

From the analysis of Fig. 6 we see that solution with 
the properties (72) is the most effective on the interval 
|cj (2) ,Cj 5) j and solution with properties (73) - on the 

interval (cj 5) , 8.oj As an example, further, the results 

of synthesis of given amplitude DP F(s l ,s 2 ) = l at 
Cj = 7.0 on the beam c 2 = 0.8q corresponding to 
different types of branching-off solutions are presented. 



The solution with the properties (73) is the most 
effective at q = 7.0 . The synthesized amplitude DP 

/(.Si,^) an d phase DP arg/(5,,5 2 ) °dd on s { , are 
given in Fig. 7, a and Fig. 7,b, respectively. Function of 
distribution of oscillating velocity in aperture (Fig. 8), 
which forms this DP, is real, but asymmetric 
concerning the axis OX . 

In Fig. 9 and Fig. 10 the results of the synthesis of a 
given amplitude DP F(s 1 ,s 2 ) = l at q = 7.0 , 
corresponding to the phase DP with properties (72), 
are presented. 

The most effective solutions of the synthesis of a given 
amplitude DP F(s 1 ,s 2 ) = l at q=7.0, corresponding 

to phase DP with properties (73), are shown in Fig. 11 
and Fig. 12. 






a) 



b) 



FIG. 7 SYNTHESIZED AMPLITUDE (A) AND ODD PHASE (B) 
DPS 

In Fig. 13 we present the results of synthesis of the 
given four-petalous DP F(s l ,s 2 ) = |sin(7ii[)| -^^(toj)! 
at q = c 2 = 7.0 , which is responsible for different 
types of solutions. The first type of solution has odd 
phase DP (Fig. 14,a) and asymmetric relatively to axis 
OY but real amplitude distribution of the oscillating 
velocity in aperture. Synthesized DP corresponding to 
this solution is shown in Fig. 13,b. 




FIG. 8 THE OPTIMAL DISTRIBUTION OF OSCILLATING 
VELOCITY, THAT FORMS, AMPLITUDE AND PHASE DP'S 
PRESENTED IN FIG. 7 
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a) b) 

FIG. 9 THE SYNTHESIZED AMPLITUDE (A) AND ODD ON BOTH 
ARGUMENTS PHASE (B) DP'S 




FIG. 10 THE OPTIMAL DISTRIBUTION OF OSCILLATING 
VELOCITY, THAT FORMS, AMPLITUDE AND PHASE DP'S 
PRESENTED IN FIG. 9 




FIG. 11 SYNTHESIZED AMPLITUDE (A) AND EVEN ON BOTH 
ARGUMENTS PHASE (B) DP'S 




FIG. 12 THE OPTIMAL DISTRIBUTION OF OSCILLATING 
VELOCITY, THAT FORMS, AMPLITUDE AND PHASE DP'S 
PRESENTED IN FIG. 11 
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a) b) 

FIG. 13 THE PRESCRIBED (A) AND SYNTHESIZED (B) 
AMPLITUDE DP'S 




a) b) 

FIG. 14 SYNTHESIZED PHASE DP (A) AND THE OPTIMAL 
DISTRIBUTION OF OSCILLATING VELOCITY IN APERTURE (B) 

The second type of solution has phase DP even on 
both arguments (Fig. 15,b) and amplitude distribution 
of the oscillating velocity in aperture symmetric 
relatively two axes (Fig. 15,a). 

Comparing the synthesized amplitude DP's (Fig. 13,b 
and Fig. 15,a) that coincide with each other to with in 
graphic precision and correspond to different types of 
solutions, we make sure that one and the same (or 
close to each other the amplitude DP) can be formed 
by different distributions of oscillating velocity in the 
aperture of antenna. 

Numerical examples of synthesis of prescribed 
amplitude funnel-shaped DP defined by formula (65) 
in domain Q are shown on Fig. 17 and Fig. 18. The 
lines of possible branching of solutions of (31) are 
given in Fig. 5. The prescribed DP and optimal 
synthesized DP are shown in Fig. 17, a and Fig. 17, b, 
respectively, at q = 9.25 , c 2 = 7.4 . The optimal 
function of distribution of oscillating velocity U (x, y) 
corresponding to synthesized DP shown in Fig. 17, b is 
presented in Fig. 18. 
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a) b) 
FIG. 15 SYNTHESIZED AMPLITUDE (A) AND PHASE (B) DP'S 




FIG. 16 THE OPTIMAL DISTRIBUTION OF OSCILLATING 
VELOCITY THAT FORMS, AMPLITUDE AND PHASE DP'S 
PRESENTED IN FIG. 15 





«) 



b) 



FIG. 17 THE PRESCRIBED (A) AND SYNTHESIZED (B) 
AMPLITUDE DP'S 




FIG. 18 THE OPTIMAL DISTRIBUTION OF OSCILLATING 
VELOCITY THAT FORMS, AMPLITUDE AND PHASE DP'S 
PRESENTED IN FIG. 17 

Conclusions 

■ Note that the presence of nonuniqueness of 
solutions in nonlinear synthesis problems of 



antennas according to the prescribed 
amplitude DP greatly complicates the 
mathematical investigations of optimal 
solutions. It provides an opportunity in 
practice to choose the most effective solution 
that has a relatively simple physical 
realization. 

■ The results of synthesis confirming the 
effectiveness of branching-off solutions are 
shown in Fig. 6. Here on the edge of segments 
5]Cj and 5 2 q functional takes the same 
values on the primary (real) and branching- 
off (complex) solutions. From here follows 
that efficiency of the synthesis on the 
branching-off solution in comparison with the 
efficiency of synthesis on the real solution can 
be achieved at smaller values of the 
parameter q characterizing the linear 

dimension of aperture. Thus, the linear 
dimension of the antenna can be decreased on 
the value Sjq or 5 2 Cj at realization of 
branching-off solution. In addition the area of 
antenna is reduced accordingly. It 
significantly influences on technical 
characteristics of antenna and the expenses of 
its realization. 

■ In the case of given amplitude DP F(s l ,s 2 ) 
even on both arguments the investigations of 
three-dimensional type of branching of 
solutions of (20), (for example, the point 
( Cl (5) ,p Cl (5) ), see Table 1), is advisable to carry 

out separately in classes of even and odd 
functions (where there is two-dimensional 
type of branching) using the invariance 
property of nonlinear integral operators 
B = (B 1 ,B 2 ) in (20) concerning the type of 

parity of phase DP. 
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